1, 2008 22:34 WSPC/INSTRUCTION FILE ccfm 



Modern Physics Letters A 

© World Scientific Publishing Company 



fct- factorization and CCFM 
the solution for describing the hadronic final states - everywhere ? 

Hannes Jung 

Department of Physics, Lund University, Box 118 
Lund, 22100 
Sweden 
hannes. jung@desy. de 

Received (Day Month Year) 
Revised (Day Month Year) 

The basic ideas of fct-factorization and CCFM parton evolution is discussed. The unin- 
tegrated gluon densities, obtained from CCFM fits to the proton structure function data 
at HERA are used to predict hadronic final state cross sections like jet production at 
HERA, but also comparisons with recent measurements of heavy quark production at 
the Tevatron are presented. Finally, the fct-factorization approach is applied to Higgs pro- 
duction at high energy hadron hadron colliders and the transverse momentum spectrum 
of Higgs production at the LHC is calculated. 
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1. Introduction 

The theory of strong interaction, QCD, has been very successful in describing many 
experimental measurements, but a number of problems have not yet been solved, 
some of which are related to the transition between the perturbative and non- 
perturbative region. The perturbative methods, however, work surprisingly well, 
even down to very low scales, where the running coupling constant, a s , starts to 
become large. Another type of problem is related to the observation, that at high 
energies, even for small values of coupling constant a s , the phase space for parton 
emissions increases fast, and that therefore it is not sufficient to include only a few 
calculable terms in the perturbative expansion. This problem can be treated by 
resumming the leading logarithmic behavior of the cross section. The most impor- 
tant contribution at small x is gluon bremsstrahlung, with the typical behavior of 
being largest in the infrared and/or collinear region. Two different resummation 
strategies have been developed: the DGLAP approach, resumming leading 

logarithms of ratios of subsequent virtualities also called collinear approach, and 
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the BFKL |SJ|H1[7] approach, resumming the infrared contributions, also called the 
^-factorization (HUH or the semi-hard approach QUEI]. The CCFM [T2[H[TI|G2] 
approach attempts to cover both the collinear and the infrared regions by con- 
sidering color coherence effects, and in the limit of asymptotic energies is almost 
equivalent [Tfilli7lll8| to the BFKL and DGLAP evolution equations. The LDC ap- 
proach [ilpUll^l^. which is a reformulation of CCFM, a unified DGLAP-BFKL 
approach 23 24 25 and the approach of doubly unintegrated parton distribu- 
tions |26j are not further discussed in this report. 

In the collinear approach the cross section is factorized into a process dependent 
hard scattering matrix element convoluted with universal parton density functions. 
Since strong ordering in virtualities is required in the evolution, the largest virtual- 
ity is in the hard scattering and therefore the virtuality of the partons entering the 
hard scattering matrix element can be neglected and they can be treated as being 
collinear with the incoming hadron a . Any physics process in the fixed order collinear 
factorization scheme is then calculated by a convolution of a process dependent co- 
efficient function C a (|) with collinear (independent of k t ) parton density functions 
at a scale [i 2 , (e.g. ^ = Q 2 in F 2 (x, Q 2 ) in deep inelastic scattering (DIS)): 



While the DGLAP approach with fixed order coefficient- and splitting functions 
is phenomenologically successful for inclusive quantities like the structure function 
F2(x, Q 2 ) in DIS, it is not fully satisfactory from a theoretical point of view, because 
"the truncation of the splitting functions at a fixed perturbative order is equivalent 
to assuming that the dominant dynamical mechanism leading to scaling violations 
is the evolution of parton cascades with strongly ordered transverse momenta" as 
Catani argued in Ref. EHI 

In the fct-factorization approach, the partons along the parton ladder are no 
longer ordered in transverse momentum. At large energies (small x) the evolution of 
parton densities proceeds over a large region in rapidity Ay ~ log(l/a;) and effects 
of finite transverse momenta of the partons may become increasingly important. 
Cross sections can then be k t - factorized [10lllll lrfll5| into an off mass-shell (kt 
dependent) partonic cross section a(^,kt) and a kt - unintegrated parton density 
function !F(z, kt): 



The unintegrated gluon density !F(z, kt) is described by the BFKL evolution equa- 
tion in the region of asymptotically large energies (small x) . 

An appropriate description valid for both, small and large x, is given by the 
CCFM evolution equation, resulting in an unintegrated gluon density A(x, kt,q), 

a However neglecting the transverse momentum of the partons even in the collinear approach has 
been criticized in Ref. 27J as being unnecessary and unphysical. 




(1) 




(2) 
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which is a function also of the additional evolution scale q described below. This 
scale q is connected to the factorization scale /i/ in the collinear approach. 

Carrying out the kt integration in eq. (J2J explicitly, a form fully consistent with 
collinear factorization can be obtained : the coefficient functions and also 

the DGLAP splitting functions leading to f a (z, lij) are no longer evaluated in fixed 
order perturbation theory but supplemented with the all-order resummation of the 
a s log 1/x contribution at small x. 

2. The CCFM evolution equation 




(a) (b) 

Fig. 1. Kinematic variables for multi-gluon emission in leptoproduction (a.) and hadroproduction 
(b). The t-channel gluon four- vectors are given by ki and the gluons emitted in the initial state 
cascade have four- vectors pi. The maximum angle (a function of the rapidity) for any emission is 
obtained from the quark box, as indicated with 3. 

The pattern of QCD initial-state radiation in a small- X event in ep and pp 
collisions is illustrated in Fig. ^ together with labels for the kinematics. According 
to the CCFM evolution equation, the emission of partons during the initial cascade 
is only allowed in an angular-ordered region of phase space. In terms of Sudakov 
variables the quark pair momentum is written as: 



p q +p q = ?(p w +EpW) + Q t 



(3) 
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where p^ and p^ are the four-vectors of incoming particles (electron-proton or 
proton-proton) , respectively and Q t is the transverse momentum of the quark pair 
in the center of mass frame of p^ and p^ (cms). The variable 3 is related to the 
rapidity Y in the cms via 



2 °\E-P Z J 2 

since E+p z = 2TyG, E-p z = 2TE^/s and E = with s = (pW +p^) 2 being 

the squared center of mass energy. Therefore S can be used to define the maximum 
allowed angle in the evolution. The momenta pi of the gluons emitted during the 
initial state cascade are given by (here treated massless): 

ft=«i(p (1) +&P (a) )+Pw, & = 4>> ( 5 ) 

svf 

with Vi = (1 — Zi)xi—i and Xi = ZiXi-\. The variables Xi and Vi are the momentum 
fractions of the exchanged and emitted gluons, while z% is the momentum fraction 
in the branching (i — 1) —* i and pu is the transverse momentum of the emitted 
gluon i. Again the rapidities yi are given by y$ = — 0.5 log ^ in the cms. 

The angular-ordered region is then specified by (Fig.^i and the lower part of the 
cascade in Fig.HJ>, for the upper part the variables have to be changed accordingly): 

£o < il < ■ ■ ■ < in < S (6) 

which becomes: 

Zi-iqi-i < Qi (7) 
where the rescaled transverse momentum qi of the emitted gluon is defined by: 

Qi = \/sii = 1 3 (8) 

The CCFM equation for the unintegrated gluon density can be written [T51I3UI 
I^Tll^j as an integral equation: 

— J ■^2 < 3>{q~zq)A s (q,zq)P gg {z,q,k t )A^- 1 k' t ,qJ (9) 

with k' t = \kt + (1 — z)q\ and q being the upper scale for any emission: 

q > z n q ni q n > z n -iq n -u Qi > Qa (10) 
The Sudakov form factor A s is given by: 

AsMo)=ex H~4 Wo dz —* ) 

with a s = ^Sa., For inclusive quantities at leading- logarithmic order the Sudakov 
form factor cancels against the 1/(1— z) collinear singularity of the splitting function. 



I 



I 
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The splitting function P gg for branching i is given by: 

p i u \ _ ®s(kti) A , , s a s (p ti ) 

igg Qij K ti ) — £-±ns \Zi,qi, rCti ) T 



(12) 



Z{ l Zr L 

with p ti = qi(l — Zi) and the non-Sudakov form factor A ns defined as: 

f 1 dz' f dq 2 

log & ns (Zi,q u k t i) = -a s J —J —<d(k u -q)<d(q-z'qi) 

The upper limit of the z' integral is constrained by the functions in ea. (|13fl by: 
Zi < z' < min(l, ku/qi), which results in the following form of the non-Sudakov 
form factor ^24s: 



(13) 



log A r 



-a s (k u )log ( — I log . 

Zi J \zoZtq, 



K ti 



(14) 



where 

1 if ku/qi > 1 
k t i/p tl if Zi < ku/qi < 1 
Zi if ku/qi < Zi 

The non-Sudakov form factor can be written as: 



1 + a s log 




+ ... + 



|ra s log(f) log(^ 



where the similarity with the Sudakov form factor becomes obvious. Note however, 
that the Sudakov form factor A s resums the large z contributions, whereas the 
non-Sudakov form factor A„ s resums the small z ones. 

In the CCFM approach the scale q (coming from the maximum angle) can be 
written as (see eq.©): 

q 2 = T 2 Es = s + Q 2 (15) 

with s = (p q + Pq) 2 and the relation of q to a particular choice of the factorization 
scale [if in the collinear approach becomes obvious. 

2.1. Improvements to the CCFM splitting function 

Originally, the CCFM splitting function P gg , as given in ea.p2|l. included only the 
soft and collinear singular terms. In the asymptotic region of large energies this 
is a reasonable choice, but effects from non-leading contributions are expected at 
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presently accessible energies. For simplicity also the scale in a s was chosen differently 
in the small and large z part of the splitting function (see en. (|12fl ^l. In Ref. 1321 it 
was suggested to use \i r = pt = q(l — z) everywhere and to include the non-singular 
terms in the splitting function. These changes are non-trivial as they need to be 
considered both in the Sudakov and non-Sudakov form factors. 

In the following improvements to the originally proposed CCFM splitting func- 
tion are discussed. 



2.1.1. The soft region 

In a DGLAP type evolution with the transverse momenta of the gluon propagators 
increasing from the proton towards the hard scattering, the non-perturbative region 
with kt < kf Ut has an influence on the initial parton density only. In CCFM, due 
to angular ordering a kind of random walk in the propagator gluon k t can be 
performed. Even during the evolution the non-perturbative region can be entered 
for k t < fe£ . In the region of small kt, ct s and the parton density are large, and 
collective phenomena, like gluon recombination or saturation might play a role. 
Thus, the fast increase of the parton density and the cross section is tamed. Much 
effort was put recently into a more detailed understanding of this special region of 
phase space (for example see Ref. [SHEIOIiDOIEl)- However, for the calculation of 
the unintegrated gluon density presented here, a simplified but practical approach 
is taken: a s (/x) is fixed for \i < Q and no emissions are allowed until k t > kf vt is 
reached. 



2.1.2. The non-singular terms in P gg 

The Pgg splitting function used in the collinear approach contains also non-singular 
terms. Such non-singular terms can be included in the CCFM splitting function, but 
care has to be taken, which terms A„ s is acting on to ensure positivity of P gg |32| : 

„ ./(1-z) z(l-z)\ A , ,f z z(l-z)\ , . 

Pgg = Ci s (k t ) + A ^T^) A - + ^) {— + 2~^j (l6) 

As the splitting function is also part of A ns and A s , they need to be modified 
accordingly The non-Sudakov form factor including the full splitting function 
is then given by: 

log A„ s = -a s (k t ) j dz + Z{l 2 Z) ^j J^-Q(kt- q)®(q - zq) (17) 



2.1.3. The scale in a s 

It was suggested in Ref. [32 to change the scale in a s to pt = q(l — z) also in the 
small z part of the splitting function P gg : 

Pgg = ^A nS + (18) 



I 



I 
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As a consequence, the non-Sudakov form factor, changes from eg . 1(14(1 to |32II37| : 



It is obvious, that a special treatment of the soft region is needed, because q can 
become very small, even q < A qcc [ at small values of z. However, in practical appli- 
cations we observe only a small effect when changing the scale of the small z part 
from k t to pt ■ 

2.2. The unintegrated gluon density 

The CCFM evolution equations have been solved numerically J50J using a Monte 
Carlo method b . The unintegrated gluon density at any x, kt and scale q is obtained 
by evolving a starting gluon distribution from the scale Qo according to CCFM to 
the scale q. The normalization N of the input distribution as well as the starting 
scale Qo, which also acts as a collinear cutoff to define z max = 1 — Qo/q, need to 
be specified. These parameters were fitted such that the structure function F 2 as 
measured at HI |39I4U| and ZEUS |41I42| can be described after convolution with the 
off-shell matrix element in the region of x < 5 • 10 -3 and Q 2 > 4.5 GeV 2 . According 
to the discussion in the last section, the following sets of CCFM unintegrated gluon 
densities are obtained: 

• JS2001 (Jung, Salam ED|) 

The splitting function P gg of ea. (|12|) is used, with Qo = 1.4 GeV. The soft region 
is defined by k c t ut = 0.25 GeV. 

• J2003 set 1 [13] 

The CCFM splitting function containing only singular terms fea. (|12H % l is used, 
with kf ut = Qo fitted to kf ut = Q = 1.33 GeV. 

• J2003 set 2 

The CCFM splitting function (ea.l|lfi|l) containing also the non singular terms is 
used. The Sudakov and non-Sudakov form factors were changed accordingly. The 
collinear cut is fitted to Q = k c t ut = 1.18 GeV. 

• J2003 set 3 gSj 

CCFM splitting function containing only singular terms but the scale in a s is 
changed from k t to pt for the 1/z term. The collinear cut is fitted to Qo = kf ut = 
1.35 GeV. The problematic region in the non-Sudakov form factor in ea. 114(1 is 
avoided by fixing a s (n) for p, < 0.9 GeV. 

A comparison of the different sets of CCFM unintegrated gluon densities is shown in 
Fig. [3 It is clearly seen, that the treatment of the soft region, defined by k t < fcf* 
influences the behavior at small x and small k t . After convolution with the off- 
shell matrix elements, all sets describe the structure function F2(^,Q 2 ) reasonably 




(19) 



b A Fortran program for the unintegrated gluon density xA{x, kt, q) can be obtained from Ref 38 
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Fig. 2. Comparison of the different sets of unintegrated gluon densities obtained from the CCFM 
evolution as described in the text. In (a — c) the unintegrated gluon density is shown as a function 

of x for different values of kt at a scale of q = 10 Ge V. In (d — f) the ratio R - 
a function of x for different values of kt is shown. 



well. In Fig. results of the fits are compared with measurements of the structure 
function ^(x, Q 2 ) as obtained by the HI 00] and ZEUS 021 collaborations. In 
Tab. H the parameters of the CCFM unintegrated gluon densities are summarized 
indicating also the quality of the fits. It is interesting to note, that the quality of 
the CCFM fits is similar to that obtained from global fits in the collinear DGLAP 
approach (e.g. Ref . llH . 



3. Comparison with hadronic final state data 

A comparison of measurements of hadronic final state properties, like jet or heavy 
quark cross sections, with theoretical predictions requires a detailed simulation of 
the experimentally accessible phase space. Such simulations are provided by Monte 
Carlo event generators, which also allow to apply the hadronization step. Monte 
Carlo event generators for DGLAP type collinear factorized processes are widely 
used (e.g. Pythia Rapgap g7j, Herwig Two Monte Carlo generators 



I 



I 
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Table 1. Parameters of the CCFM unintegrated gluon densities. The x / n df 
is obtained from a comparison to HERA F2 data 39 40 41 42 in the range 
as < 5 • lCr 3 and Q 2 > 4.5 GeV 2 . 



set 






Qo (GeV) 


kf ut (GeV) 


X 2 /ndf 


JS2001 3D 


ca.rm 


ca.tm 


1.40 


0.25 


1197/248 = 4.8 


J2003 set 1 


cq.JEJ 


ea.lul 


1.33 


1.33 


321/248 = 1.29 


J2003 set 2 


eg. 1161 


eq.JTTJ 


1.18 


1.18 


293/248 = 1.18 


J2003 set 3 


ca.llSl 


eq.flU 


1.35 


1.35 


455/248 = 1.83 



have been developed including the small x evolution equations: Cascade |3UII38| . 
which follows explicitely the CCFM approach of angular ordering and Ldc |19l2UI2l1 
122] which is a reformulation of CCFM (not described here). It has been shown |30j . 
that the parton shower approach used in Cascade reproduces exactly the properties 
of the CCFM evolution described in section |2J The Cascade Monte Carlo event 
generator has been frequently used for comparison with HERA measurements, like 
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heavy quark [49 , 50 , 51 h% and high p t jet [F3 , 30 , 54 55 production, and also with 
bottom production at the Tevatron |52j . 

In the following sections a few examples are presented where data arc com- 
pared to predictions obtained with Cascade based on the new unintegrated gluon 
densities, described in section 12.21 

3.1. Jet cross section at HERA 

The azimuthal correlation of dijets at HERA is sensitive to the transverse momen- 
tum of the partons incoming to the hard scattering process and therefore sensitive 
to the details of the unintegrated gluon density. This was studied in a measure- 
ment [S3] of the cross section for dijet production with Et > 5(7) GeV in the 
range 1 < rji a {, < 0.5 in deep-inelastic scattering (10 -4 < x < 10~ 2 , 5 < Q 2 < 100 
GeV 2 ). In LO collinear factorization, dijets at small xsj are produced essentially 
by 73 - * 99, with the gluon collinear to the incoming proton. Therefore the qq 
pair is produced back-to-back in the plane transverse to the j*p direction. From 
NLO (0(a 2 )) on, significant deviations from the back- to-back scenario can be ex- 
pected. In the fc t -factorization approach the transverse momentum of the incoming 
gluon, described by the unintegrated gluon density, is taken explicitely into account, 
resulting in deviations from a pure back-to-back configuration. The azimuthal de- 
correlation, as suggested in Ref. ESI can be measured: 

S = J ° 5 f_j_^_J <a<180 (20) 

r 180 iv 2 _, et (A0*,x,g 2 )dA0* 

In the measurement shown in Fig. [3J a = 120° has been chosen. The data are 
compared to predictions from Cascade using J 2003 set 1-3. Also shown for 
comparison is the NLO-dijet \57\ calculation of the collinear approach. One clearly 
sees, that a fixed order NLO-dijet calculation is not sufficient, whereas J2003 set 
2 gives a good description of the data. However, the variable S is sensitive to the 
details of the unintegrated gluon distribution, as can be seen from the comparison 
with J2003 set 1 and set 3. 

A measurement, aiming to observe deviations from the collinear DGLAP ap- 
proach, is the production of jets in the forward (proton) region. The phase space is 
restricted to a region of Q 2 > 5 GeV 2 and E 3 T et > 3.5 GeV in the forward region 
of 1.7 < rjj e t < 2.8 with the additional requirement of 0.5 < Ej,- et /Q 2 < 2, a re- 
gion where the contribution from the evolution in Q 2 is small. The cross section for 
forward jet production has been measured by HI |5.'i| as a function of xbj, shown 
in Fig. [SI together with predictions from Cascade using J2003 set 1 - 3. Also the 
NLO-dijet [37J prediction in the collinear approach is shown. The fixed NLO-dijet 
calculation falls below the measurement, whereas the fct-factorization approach sup- 
plemented with CCFM evolution gives a reasonable description of the data for all 
J2003 set 1 - 3. 



1, 2008 22:34 WSPC/INSTRUCTION FILE ccfm 



11 




X 


250 


T3 




6 225 


T3 






200 




175 




150 




125 




100 




75 




50 




25 








• H1 preliminary 

jets p, > 3.5 GeV 



J2003 - set 1 
J2003 - set 2 
J2003 - set 3 
NLO dijet 



0.001 



0.002 



0.003 



0.004 



Fig. 4. The S distribution for dijet events 
1541 . The prediction from CASCADE using 
J2003 set 1 and J2003 set 2 are shown to- 
gether with the NLO-dijet prediction in the 
collincar factorization approach. 



Fig. 5. The cross section for forward jet 
production as measured by HI 1531 as a func- 
tion of xgj . The prediction are the same as 
in Fig.|3] 



3.2. Charm Production at the Tevatron 

The differential cross section as a function of the transverse momentum of D-mesons 
has been measured in pp collisions at y/s = 1.96 TeV by the CDF collaboration 
They find the measured cross section to be larger than the NLO predictions in the 
collinear factorization approach by about 100 % at low pt and 50 % at high p t . 
In Fig. El the measurement is shown together with the predictions obtained in fc t - 
factorization using Cascade with the CCFM unintegrated gluon density described 
above. For J2003 set 1 and set 3 good agreement for all measured charmed mesons is 
observed. The cross section predicted using J2003 set 2 falls below the measurement, 
in contrast to the observation in Fig. showing the sensitivity of the measurements 
on the details of the unintegrated gluon density. 



3.3. Bottom Production at the Tevatron 

The cross section for bb production in pp collision at y/s = 1800 GeV has also 
been compared with the prediction of Cascade based on the CCFM gluon den- 
sities. Since Cascade generates full hadron level events, direct comparisons with 
measured cross sections for the production of b quarks decaying semi-leptonically 
into muons are possible. The muon cross sections as a function of the transverse 
momentum p^ and pseudo-rapidity \y p, \ as measured by DO |59| are compared to 
the Cascade prediction in Fig. [7| and Fig. El Both the and \y^\ cross sections 
are well described with J2003 set 1 and set 3. The prediction based on J2003 set 
2 falls below the measurement, as was already observed in the charm case, again 
indicating the sensitivity to the details of the unintegrated gluon density. 
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Fig. 6. Differential cross section for D meson production as measured by CDF 1581 as a function 
of the transverse momentum compared to predictions from CASCADE. 



3.4. Higgs production in pp 

Higgs production at high energies proceeds predominantly via gluon - gluon scat- 
tering. Much effort has been put into the calculation of higher order corrections in 
the collinear approach, not only for the calculation of the total production cross 
section, but also for the calculation of the transverse momentum spectrum of the 
Higgs boson [6()lf>ll(S2] . At the large Tevatron or LHC energies and with an expected 
Higgs mass of O(100 — 200) GeV, the £; t -factorization approach can be also used 
to estimate higher order corrections. The off-mass-shell matrix element g*g* — > h 
has been calculated in Ref. 63 in the high energy approximation for m t — > oo and 
is now implemented in Cascade. In Fig. the distribution of the longitudinal 
momentum fractions x of the gluons and their transverse momenta are shown for 
Tevatron and LHC energies. In both cases, the longitudinal momenta reach values 
of the same order as the average transverse momenta (~ 10(18) GeV for Tevatron 
(LHC) energies, respectively), making the fc t -factorization approach applicable. It 
can be seen in Fig. [5J that the kt spectrum of the gluons is different for the different 
unintegrated gluon densities described in sec. 12.21 both in the small and also large 
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Fig. 7. Cross section for muons from b- 
quark decays as a function of p/^, (per unit 
rapidity) as measured by DO 1591 compared 
to the prediction of CASCADE. 
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Fig. 8. Cross section for muons from b- 
quark decays as a function of for two 
different p^, cuts as measured by DO 1591 
compared to the prediction of CASCADE. 



Table 2. Cross section for Higgs boson 
production cr(pp — > hX) 



gluon distribution 








1.96 TeV 


14.0 TeV 


J2003 set 1 


0.47 pb 


22.6 pb 


J2003 set 2 


0.38 pb 


7.8 pb 


J2003 set 3 


0.54 pb 


24.6 pb 



kt region. In Fig. 1101 the differential cross section for Higgs production as a function 
of pt obtained with Cascade is shown for the different unintegrated gluon densi- 
ties. The total cross section for Higgs production is given in Tab. [21 It is interesting 
to note, that the different set of unintegrated gluon densities predict similar cross 
sections as a function of pt for Tevatron energies whereas at LHC the cross section 
differs by factors up to 3. This again clearly indicates the sensitivity to the details 
of the unintegrated gluon density, which can be determined much more precisely 
with the forthcoming measurements at HERA. A similar result has been obtained in 
Ref. |M] 1651 approximating the matrix element and the CCFM unintegrated gluon 
distribution and also applying doubly unintegrated parton distributions in Ref. 



4. Conclusion 

It has been shown, that /c t -factorization and the CCFM evolution of the gluon den- 
sity is a powerful tool for the description of hadronic final state measurements. The 
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Fig. 9. Longitudinal and transverse momenta of the gluons for g* g* — > h for Tcvatron 
(upper) and LHC (lower) energies. 



unintegratcd gluon densities, which are obtained from CCFM evolution convoluted 
with the off-mass-shell matrix elements to describe HERA F 2 data are implemented 
in the full hadron level Monte Carlo event generator Cascade. 

Jet measurements at HERA, but also measurements of charm and bottom pro- 
duction at the Tevatron can be reasonably well described (with specific sets of 
unintegrated gluon densities), whereas calculations performed in the collinear ap- 
proach even in NLO have difficulties to describe the data. This shows the advantage 
of applying /c t -factorization to estimate higher order contributions to the cross sec- 
tion but also the importance of a detailed understanding of the parton evolution 
process. 

The same unintegrated gluon densities have been used to calculate also the trans- 
verse momentum spectrum of the Higgs boson production at the Tevatron and the 
LHC. Given the large transverse momenta of the gluons involved, ^-factorization is 
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Fig. 10. Differential cross section for Higgs production as a function of the transverse momentum 
for muiggs = 125 GeV obtained with the different un-integrated gluon densities. 

the appropriate tool for calculating higher order corrections. However, the different 
unintegrated gluon densities show significant effects at LHC energies, indicating the 
need for better experimental constraints as well as further theoretical studies for a 
more detailed understanding of parton evolution at large energies. 
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